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1.   Introduction. 

The  bivariate  delayed  Poisson  process  is  a  particular  type  of  bivariate 
point  process,  that  is  a  point  process  with  two  types  of  event  occurring  along 
a  time  axis.   Bivariate  point  processes  were  first  studied  by  Cox  and  Lewis 
(1972)  in  their  Berkeley  Symposium  paper  and  here  their  analysis  of  the 
bivariate  Poisson  processes  (without  Poisson  noise)  is  carried  further.   This 
bivariate  point  process  is  formed  by  subjecting  the  events  of  a  main  Poisson 
process  to  independent  pairs  of  delays;  the  two  types  of  delay  for  each  main 
event  produce  the  two  types  of  event  in  the  bivariate  point  process.   Since  it 
is  well  known  that  a  Poisson  process  whose  events  are  independently  and  identi- 
cally displaced  remains  a  Poisson  process,  the  events  of  each  type  considered 
separately  form  two  Poisson  processes.   They  cannot,  however,  be  independent 
Poisson  processes  because  their  events  are  associated  in  pairs  by  the  displace- 
ment centres  (the  events  of  the  original  Poisson  process).   The  resulting 
bivariate  point  process  is  by  definition  a  particular  type  of  bivariate 
Poisson  process  because  its  marginal  processes  for  events  of  each  type  are 
themselves  (univariate)  Poisson  processes.   Cox  and  Lewis  (1972)  assumed  a 
model  in  which  the  delays  within  pairs  were  independent  and  non-negative,  and 
discussed  the  transient  and  asynchronous  bivariate  counting  processes.   We  shall 
mainly  assume  that  the  delays  within  pairs  are  dependent  and  two-sided  and  will 
first  obtain  the  stationaly  initial  conditions.   These  lead  to  the  joint  distri- 
bution of  the  numbers  of  delays  of  each  type  in  progress  at  an  arbitrary  time 
and  a  general  form  of  the  asynchronous  bivariate  counting  distribution.   Finally, 
we  obtain  the  semi-synchronous  marginal  interval  distribution,  that  is  the 
distribution  of  the  time  from  an  event  of  one  type  to  the  next  event  of  the 
opposite  type;  such  distributions  are  of  obvious  use  in  the  statistical  analysis 


of  bivariate  point  processes,  and  have  been  studied  in  the  context  of  the 
general  theory  by  Wisniewski  (1972).   Milne  (unpublished  Ph.D.  thesis),  has 
shown  that  the  bivariate  Poisson  process  studied  here  is  the  most  general  form 
for  an  infinitely  divisible  bivariate  point  process. 

2.  Notation  and  Description  of  the  Model. 

The  bivariate  delayed  Poisson  process  is  constructed  from  a  main  Poisson 

process  of  rate  u.  Associated  with  each  of  its  events  is  a  pair  of  delays 

Y   and  Y  ;   thus  an  event  at  time  t  in  the  main  Poisson  process  produces 

an  event  at  time   t  +  Y   in  the  a-process  and  an  event  at  time  t  +  Y,   in 

a         — c b 

the  b-process.   We  will  assume  that  the  random  variables  Y   and  Y,   with 
— c a       b 

ranges   (-00,00)  have  finite  means  and  a  dependent  joint  distribution  without 
any  atoms.   The  pairs   (Y  ,Y  )  associated  with  different  main  events  are 
assumed  to  be  independently  and  identically  distributed. 

The  random  variables  N  (t,t+t  ),   N,  (t,t+t,  )  denote  the  numbers  of 
type-a  and  type-b  events  in  the  intervals   (t,t+t  ]   and   (t,t+t,  ],   respectively; 

3.  D 

when  there  can  be  no  confusion  these  will  be  abbreviated  to  N  (t  )  and  N,  (t,  ), 

a  a        b  b 

with  joint  probability  generating  function   (j  .p.g.f  .),(p(z- ,z„;  t  , t,  ). 

3.  Stationary  Initial  Conditions. 

The  purpose  and  methodology  of  stationary  initial  conditions  for  uni- 
variate point  processes  have  been  described  in  Lawrance  (1972);  they  carry 
over  quite  obviously  to  bivariate  processes  and  so  will  not  be  laboured  here. 
We  wish  to  study  the  behaviour  of  the  process  at  an  arbitrary  time,  and  this 
is  defined  as  the  time   t  as   t  -*■  °°.   It  is  apparent  for  the  bivariate  delayed 
Poisson  process  that  at  any  time  t   there  may  be  a  (random)  number  of  events 
of  each  type  in  the  course  of  being  delayed.  The  initial  conditions  must 


therefore  specify  jointly  the  distribution  of  these  numbers  and  the  individual 
residual  delay  times  measured  from  t.   We  begin  by  dealing  with  the  numbers 
and  delay  times  which  are  associated  with  main  events  occurring  before  time 
t;  however,  when  delays  can  be  negative  (or  advances),  we  have  also  to  consider 
delays  associated  with  main  events  after  time  t. 

The  behaviour  at  time  t  yields  to  an  analysis  of  the  pairwise  behav- 
iour of   (Y  ,Y,  )   for  each  main  event  in   (0,t].   There  are  four  possibilities 
a  b 

concerning  the  type-a  and  type-b  events  which  terminate  Y   and  Y  :   both 
occur  after  time   t,   type-a  occurs  after  time   t  and  type-b  before  time   t, 
vice-versa,  and  both  occur  before  time  t.   Let  there  be  n  main  events  in 
(0,t],  with  Z  ,Z^   ,Z^b)   and  n  -  Z  -  Z^a)  -  Z^b)  main  events  giving  pairs 
of  type-a  and  type-b  events  of  the  above  four  types,  respectively.   Let 


{Y  (i),Y,  (i)}   for  i  =  1,2,...,Z   be  the  delay  times  past   t  associated 
a     b  t 


with  main  events  counted  by   Z  ,   and  similarly  let  Y'(i)   for  i  =  1,2,...,Z 

and  Y'(i)   for  i  =  1,2,...,Z     be  the  delay  times  past  t  asociated  with 
b  t 

the  other  two  situations. 

The  key  feature  of  the  present  model  which  makes  it  tractable  is  that, 
given  the  n  main  events  in  (0,t],  they  are  independent  and  uniformly  dis- 
tributed over  this  interval.   Thus  we  can  see  that 


(a) 
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In  (3.1)  the  variable  v  is  used  to  denote  that  the  typical  main  event  being 
considered  is  located  at  time  t  -  v,   or  a  time  point  v  before  t.   Letting 
t  -*•  »  in  (3.1)  we  obtain  after  a  little  simplification,  and  dropping  the 
suffix  t,   the  arbitrary  time  initial  conditions  as 


Z  =  k,Ya(i)  k  a±,  Yb(i)  is  b±  for  i  =  l,2,...,k 
Z(a)  =  k  ,  Y'(i)  k  a!   for  i  =  1,2 k 

cl      a.  X  a 

_z(b)  =  1^,  Yb(i)  *  b»   for  i-1,2,...,!^ 


=  exp{-  y 


[1-P(Y  <v,  Y,<v)]dv} 
a     b 


v=0 
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a    l   d    i 
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P(Y  >v+a»   Y,<v)dv] 
a    i   b 


v=0 


-1  K 
x  (k»)  X  Trb[y 

^     i-1 


P(Y  <v,  Y.>v+b!)dv]. 
a     b    i 


v=0 


(3.2) 


Both  (3.1)  and  (3.2)  are  valid  for  k,k  ,le  £  0  subject  to  the  interpretation 
that  when  any  of  k,k  ,le   are  zero  the  corresponding  products  are  omitted. 

We  next  consider  the  behaviour  at  any  particular  time,  which  we  take  as 
zero,  which  is  due  to  negative  delays.   We  must  consider  the  behaviour  at  0 
due  to  main  events  in   (0,t]   and  then  let   t  -*■  °°.   The  situation  is  similar 
to  that  just  considered  except  that  we  are  interested  in  the  'overshoots'  of 
the  delays  past  zero  in  the  negative  direction.   For  each  event  in   (0,t]   we 
either  have  a  double  overshoot,  an  overshoot  of  type-a  only,  an  overshoot  of 
type-b  only,  or  no  overshoot  at  all.   Let  these  overshoots,  measured  as  posi- 
tive quantities,  be  denoted  by  X  (i),  X,(i)   for  i  =  1,2, ...,k,   X' (i)   for 

3-D  Si 

i  =  l,2,...,k   and  ^/(i)   for  i  =  l,2,...,k,  .   Here  we  have  changed  from 
Y  to  X  in  our  notation  of  the  'forward'  initial  conditions,  and  will  use 
C  instead  of  Z  for  counting  the  main  events.   The  'backward'  initial  condi- 
tions can  now  be  obtained  analogously  to  the  forward  ones,  and  yield  after  some 
reduction 


k,  X  (i)  £  a  ,  X  (i)  ^  b±  for  i  =  1,2,...,] 


£(a)  =  k  ,  X'(i)  £  a!   for  i  =  l,2,...,k 
a   a       i  'a 
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=  exp{-  y 
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[1-P(Y  >v,  Y  >v)]dv} 
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-1 
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a     b    l 


(3.3) 


V=-c 


A  point  to  notice  about  both  the  forward  and  backward  initial  conditions  is 
that  they  do  not  simplify  appreciably  when  the  delays  within  pairs  are  inde- 
pendent.  Even  if  Y   and  Y   are  independent  in  the  model,  the  initial 
pairwise  delays,   X  (i)  and  X  (i),   if  there  are  any,  will  be  dependent. 

4.   The  Numbers  of  Events  in  the  Course  of  Being  Delayed. 

One  aspect  of  the  bivariate  delayed  Poisson  process,  considered  at  an 

arbitrary  time   t,   is  the  number  of  events  in  the  main  process  generated 

before  time   t  which  are  delayed  until  after  time  t,   and  conversely  the 

number  of  events  generated  after  time   t  which  are  advanced  to  before  time 

t.   Now  the  bivariate  delayed  Poisson  process  can  be  regarded  as  a  simple 

bivariate  case  of  Lewis's  branching  Poisson  process  (Lewis,  1964),  where  a 

main  event  generates  two  subsidiary  processes  each  of  one  event.   In  this 

terminology  we  will  be  obtaining  the  joint  distribution  of  the  numbers  of  each 

type  of  subsidiary  process  running  at  an  arbitrary  time.   It  is  shown  in  Lewis 

(1964)  for  positive  independent  delays  that  marginally  the  distribution  will 

be  Poisson  with  means   E(Y  ),   E(Y,  );   and  also,  for  positive  independent 

a       b 

delays,  the  joint  Poisson  distribution  is  stated  in  Cox  and  Lewis  (1972). 
The  result  of  two-sided  dependent  delays  is  derived  here,  and  perhaps  places 
the  earlier  results  in  a  wider  perspective. 

The  required  joint  distribution  is  obtained  from  the  forward  and  back- 
ward initial  conditions  by  first  setting  the  a., a!   and  b.,b!   in  (3.2)  and 

J  °  11        11 

(3.3)  equal  zero.   The  backward  and  forward  contributions  will  be  independent 
because  of  the  main  Poisson  process,  and  we  require  the  joint  distribution  of 
{Z+Z(   ,  Z+Z^bh  +  {£.+?,,£+£}.      Let  the  probability  generating  functions 
(p.g.f.'s)  of  the  two  pairs  be  cp  (z  ,z„)   and  cp  (z  ,z_).   From  (3.2)  we  thus 
have 


OO  00  00 
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7=0 


This  is  the  p.g.f.  of  a  bivariate  Poisson  distribution  X  and  is  essentially 
the  result  stated  by  Cox  and  Lewis  (1972,  equation  (4.34)).  For  the  negative 
delays    (or  advances)  we  have   from   (3.3) 


(J)2^Z1,Z2^   =  exP^(z1"1) 


0 

P(Y   <v)dv  +  y(z_-l)  P(Y<v)dv 

a  z  b 


V=— oo 
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+  y(Zl-l)(z2-l) 


V=— oo 


P(Y   <v,    Y^<v)dv}    . 
a  b 


(4.2) 


V=— oo 


The  product  of  (4.1)  and  (4.2)  gives  cp(z  ,z_),  the  required  p.g.f.,  and 
remains  of  the  well  known  infinitely  divisible  bivariate  Poisson  form.   The 
sum  of  the  first  integrals  in  (4.1)  and  (4.2)  shows  that  the  marginal  means 
are  E{ |Y  |}  and  E{ |Y  |},   respectively.   We  now  show  that  the  covariance 
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v=— °° 
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is,  apart  from  the  factor  y,   the  expectation  of  M(Y  ,Y,  )  where 

a  b 


0 


for  Y  Yv  £  0 
a  b 


M(Ya,Yb)  = 


Min{  Y  L  IyJ}   for  Y  Y  >  0, 
a     b  a  b 


(4.4) 


This  may  be  seen  from  the  successive  equalities 


E{M(Y  ,Y,  )}  =  E{min(Y  ,Y,  )  |Y  >0,  Y  >0}P(Y  >0,Y  >0) 
a  b  a  b  '  a     b      a    b 
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(4.5) 


The  interpretation  of  M(Y  ,Y,  )   is  that  it  is  the  absolute  length  of  the  over- 

a  b 

lap  of  a  pair   (Y  ,Y  ) ;   if  Y   and  Y   are  of  opposite  sign  there  is  no 
overlap  and  M(Y  ,Y  )   is  zero.   The  joint  distribution  of  the  numbers  of 
events  which  are  in  the  course  of  being  delayed  at  an  arbitrary  time  thus  has 
joint  p.g.f. 

(p(Zl,z2)  =  exp{(Zl-l)yE(|Ya|)  +  (z2-l)yE(  |Yfe  |  )  +  (Z]-l)  (z2-l)uE[M(Ya,Yb)  ]  }.    (4.6) 

When  delays  are  non-negative,   |Y   =  Y  ,  Iy^   =  Yv  ,  M(Y  Yu)  =  min(Y  ,Y,  )   and 

1  a '    a  '  b '    b     a  b        a  b 

(4.6)   reduces  the  result  stated  by  Cox  and  Lewis  (1972).   We  may  also  note 


8 


that  (4.6)  applies  to  a  bivariate  generalized  branching  Poisson  process, 

Lewis  (1969),  if  Y   and  Y,   are  understood  as  the  durations  of  subsidiary 

a        b 

processes.   Further  to  Cox  and  Lewis's  remarks  in  section  4.6  of  their 
Berkeley  paper,  it  is  interesting  to  note  that  the  initial  conditions  cannot 
be  constructed  from  the  joint  probabilities  in  (4.6)  and  the  marginal  forward 
recurrence  time  distributions  of  Y   and  Y  ;   this  is  because  the  dependency 
structure  requires  the  joint  distribution  of  these  quantities. 

5.   The  Stationary  Bivariate  Counting  Process. 

In  this  section  we  obtain  the  joint  p.g.f.  of  N  (t  )   and  R(t,  ), 
the  number  of  events  of  type-a  and  type-b  which  occur  in  the  intervals   (0,t  ] 

a. 

and   (0,t,  ]   of  the  stationary  process.   We  know  marginally  that  N  (t  )   and 
II  (t,  )   are  Poisson  random  variables  of  means   ut   and  ut,  ,   so  interest 
concerns  their  joint  behaviour.   Cox  and  Lewis  (1972,  Sections  4.3  and  4.4) 
obtained  the  joint  distribution  of  N  (t  )   and  K  (t,  )   for  pairwise  independ- 
ent and  non-negative  delays,  and  we  present  here  the  corresponding  result  for 
pairwise  dependent  and  two-sided  delays;  this  more  general  case  is  shown  to 
produce  a  rather  simpler  result. 

The  type-a  and  type-b  events  in   (0,t]   arise  from  three  sources:   (1) 
those  from  main  events  in   (-°°,0]   which  have  been  delayed  into   (0,t  ]   or 

cL 

(0,t,  ];  (2)  those  from  main  events  in   (0,max(t  ,  t,  )]  which  have  been  delayed 
b  a   b 

or  advanced  and  are  still  in   (0,t  ]   or   (0,t,  ];  (3)  those  from  main  events 

a  D 

in   (max(t  ,  t,  ),°°]   which  have  been  advanced  into   (0,t  ]   or   (0,t,  ].   The 
contributions  from  these  three  sources  are  independently  distributed  in  virtue 
of  the  Poisson  main  process,  and  we  shall  let  their  bivariate  p.g.f's  be  denoted 
by  (p(z1,z2;-»,0),   cp(z  ,z2;0,ta,tb)  and  <p(z1,  z£;  t&t  tb>») ,   respectively. 


Derivations  of  all  three  are  similar;  particular  cases  of  the  first  two  appear 
in  Cox  and  Lewis  (1972),  and  so  we  limit  ourselves  to  the  derivation  of 
(p(z-,z_;t  ,  t,  ,°°)  which  arises  from  negative  delays.  We  first  state  for 
future  reference, 


(p(z1,z2;-<»,0)   =  expfuCz^l) 


P(v<Ya<v+ta)dv  +  y(z2-l)        P(v<Yb<v+tb)dv 
v=0  v=0 


+  u(Zl-l)(z2-l) 


P(v<Y  <v+t  ,  v<Y  <v+t,  )dv>     (5.1) 
a    a     b    b 


v=0 


and 


(p(z1,z2;0,ta,tb)  =  exp{y(z1-l) 


max(ta,tb) 


P(-v<Y  <t  -v)dv  +  y(z_-l) 
a  a  2. 


v=0 


mfc(ta'tb) 

P(-v<Y  <t,-v)dv 

v=0 


max(ta,tb) 


+  y(Zl-l)(z2-l) 


P(-v<Y  <t  -v,  -v<Y  <t, -v)dv>.   (5.2) 
a  a       b  b 


v=0 


For  cp(z  ,z~;t  ,t,  ,°°)   consider  events  which  are  advanced  from   (max(t  ,t,  ),T] 
into   (0,t  ]  or   (0,t_  ],   and  ultimately  let  T  -»■  °°.   For  a  main  event  at 
time  max(t  ,t  )  +  v,   its  contribution  to  (p(z-,z_;t  ,  t,  ,°°)   is 

1  +  (z1-l)P(-v-t()<Ya<-v-t1)  +  (z2-l)P(-v-t0<Yb<-v-t2) 

+  (z1-l)(z2-l)P(-v-t0<Ya<-v-t1,  -v-tQ<Yb<-v-t2)     (5.3) 

where   t»  =  max(t  ,t,),  tn  =  max(t,-t  ,0),  t_  =  max(t  -t,  ,0). 
0       a  b    1       b  a      2       a  d 

Conditional  on  the  number  of  events  in   (tn,T],   the  terms  (5.3)  are  independ- 
ently and  uniformly  distributed  over   (t_,T].   Hence,  following  the  argument 
in  section  3  to  obtain  the  initial  conditions,  we  obtain 
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<p(z1,z2;ta,tb,»)  =  exp{y(Z;L-l)   P(-v-tQ<Ya<-v-t1)dv  +  y(z2~l)   P(-v-tQ<Yb<-v-t2)dv 

v=0  v=0 

00 

+  Vi(Z:L-l)(z2-l)    P(-v-t0<Ya<-v-t1,-v-t()<Yb<-v-t2)dv}  .        (5.4) 
v=0 

We  observe  that  (5.1),  (5.2)  and  (5.4)  are  all  of  the  infinitely  divisible 
bivariate  Poisson  form,  and  thus  it  remains  only  to  obtain  expressions  for  the 
parameters  in  their  product;  this  is  the  required  bivariate  p.g.f.   Simple 
changes  of  variables  in  the  integral  coefficients  of   (z  -1)   show  that  their 
sum  may  be  written 

+oo 

y    P(v<Y  <v+t  )dv.  (5.5) 

J       a    a 

V=— oo 

If  we  assume  that  Y   has  a  density  then  this  integral  evaluates  to  yt 

3.  3. 

straight  away;  generally  we  will  make  this  assumption.   However,  if  Y   is 

3. 

restricted  to  non-negative  values  this  may  easily  be  proved  without  the  assump- 
tion of  a  density.   We  thus  have  that  N  (t  )   and  N,  (t,  )   have  marginal 

a  a        b  b 

Poisson  distributions  of  means  yt   and  yt,  .   Similarly,  adding  together  the 
coefficients  of   (z..-l)  (z.-l)   from  (5.1),  (5.2)  and  (5.4)  gives 


+oo 

M 

V=— oo 


P(v<Y  <v+t  ,  v<Y  <v+t.  )dv  (5.6) 

a    a     b    d 


as   Cov{N  (t  ),  N,  (t,  )},   a  quantity  which  is  always  non-negative. 
Thus,  the  general  result  for  the  bivariate  p.g.f.  of  N  (t  )   and  N^CO   is 
cp(z1,z2;ta,tb)  =  exp{(z1~l)yta  +  (z^Dyt^ 

+00 

+  (z1-l)(z2-l)y    P(v<Ya<v+ta,  v<Yb<v+tb)dv}  .    (5.7) 


V=— oo 
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When  Y   and  Y   are  non-negative,  the  covariance  becomes 

a  d 


mirn(ta»tb) 


v=0 


0<Y  <t  -v 
a  a 

0<Y,  <t, -v, 

b   b  -* 


dv  +  v        P 
v=0 


v<Y  <v+t 
a    a 

Lv<Yb<v+tbJ 


dv, 


(5.8) 


and  this  makes  (5.7)  correspond  to  Cox  and  Lewis  (1972,  equation  (4.25)). 
Setting  z  =  z_  =  0  in  (5.7)  gives  the  bivariate  exponential  forward 
recurrence  times;  further  details  are  given  by  Cox  and  Lewis. 


6.   The  Semi-Synchronous  Interval  Distributions. 

Finally,  in  this  paper  we  consider  the  distributions  of  intervals 
between  events  in  the  bivariate  delayed  Poisson  process;  the  intervals  between 
events  of  the  same  type  have  trivially,  exponential  distributions  of  parameter 
u.   Consequently,  we  shall  be  interested  in  the  distribution  of  intervals 
between  consecutive  events  of  opposite  type;  these  have  been  called  by  Cox 
and  Lewis  (1972)  the  semi-synchronous  intervals.   We  take  an  arbitrary  event 
of  type-a,  say,  and  require  the  distribution  of  the  interval  to  the  first 
subsequent  event  of  type-b;  this  interval  will  be  denoted  by  X   .   Some 
general  theory  relating  synchronous,  semi-synchronous  and  asynchronous  interval 
distributions  is  given  in  Wisniewski  (1972).   We  shall,  however,  develop  our 
result  from  the  definition  in  Cox  and  Lewis  (1972), 


P(X(b)>x)  =  lim  P{R(t,t+x)  =  0  IN  (0,t)  ^  1}  , 
t+0    * 


(6.1) 


which  is  analogous  to  the  definition  in  Khintchine  (1960,  Section  9)  of  Palm 
functions  for  synchronous  interval  distributions.   Taking  (6.1)  we  write  it  as 


P(X(b)>x)  =  lim 
a       t+0 


P{KL(t,t+x)  =  0}  -  P{N  (0,t)  =  0,  ^(tjt+x)  =  0} 


P{N  (0,t)  ^  1} 

3 


(6.2) 
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Dividing  both  the  numerator  and  denominator  by   t,   and  assuming  that  the 
individual  limits  exist,  as  will  be  evident  in  a  moment  for  the  present 
process,  we  have  the  general  result 


P(xfb)>x)  =  -u"Vp{N  (0,t)  =  0,  N,  (t,t+x)  =  01  n. 
a  t   a  b  \t=(J; 


(6.3) 


Here  D   denotes  a  right  hand  derivative  and  y   is  the  value  of 

t~  P{N  (0,t)  £  1}   as   t  +  0;   the  latter  is  trivially  true  from  Khintchine's 

work  since  N  (0,t)   counts  events  in  a  Poisson  process  of  rate  y. 

3. 

Our  main  concern  now  is  thus  to  evaluate  P{ (N  (0,t)  =  0,  N,  (t,t+x)  =  0) 
for  the  bivariate  delayed  Poisson  process.   Since  N  (0,t)   and  N^(t,t+x) 
refer  to  the  stationary  process  we  must  consider  main  events  which  have  occurred 
in  the  intervals   (-T  ,0],  (0,t]   and   (t,t+T-],   eventually  letting  T   and 
T  -*■  0°.   A  main  event  which  occurs  at  -v  in   (-T  ,0]   does  not  lead  to  an  a- 
event  in   (0,t]   or  a  b-event  in   (t,t+x]  when  its  delays  are  such  that 
Y   €  [v,t+v]',   Y   €  [ t+v , t+v+x ] ' ,   where  prime  denotes  the  complementary 

3.  D 

interval;  similarly  for  a  main  event  at  v  in   (0,t]   the  condition  on  the 
delays  is  Y   €  [-v,t-v]',  Y   €  [t-v, t-v+x] ' ,   and  finally  for  an  event  at 

3  D 


t+v  in   (t,t+T  ]   the  condition  is  Y   €  [-(t+v) ,-v] ' ,  Y   €  [-v,x-v]'. 
Use  of  the  conditional  property  of  the  main  Poisson  process  for  each  of  the 


intervals   (-T  ,0],  (0,t]   and   (t.t+Tj,   leads  to  the  result 


P{N  (0,t)  =  0,  N  (t,t+x)  =  0}  =  exp{y 

3.  D 


v=0 


(Y  ^[v,t+v]'   > 
[-1WY"«[«v,t4«]'J]*' 


c 

t 


[-1+P 


v=0 


Y  €[-v,t-v]' 
a 

Yv€[ t-v, t-v+x] ' 
b 


+oo 


]dv  +  y 


[-1+P 


v=0 


rY  €[-(t+v),-v]'>, 
a 

Y,  €[-v,x-v]' 

D 


]dv}  .     (6.4) 
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Slight  changes  of  variable  in  the  second  and  third  terms  then  give 


P{Na(0,t)  =  0,  ^(t.t+x)  =  0}  =  exp{y     [-1+P 


■Ya€[u,t+u]'    > 

Y  .([t+u,t+u+x]; 


]du}  .     (6.5) 


When  the  delays  Y   and  Y   are  non-negative,  the  corresponding  result  is 

3.  D 


P{N  (0,t)  =  0,  NL(t,t+x)  =  0}  =  exp{y 
a  b 


[-1+P 


+  y 


[-1+P 


u=0 


u=0 

rY  >u 
a 

.Yb€[u,u+x]'J 


(Y  €[u,t+u]f 
Yb€[t+u,t+u+x]'J 


]du 


]du  -  y 


P(Y  <u)du)    . 
b 


(6.6) 


u=0 


We  can  now  use    (6.3)    for   the   serai-synchronous   interval   distribution  and  have 


-H» 


P{X(b)>x}   = 


u= 


[1-P 


,Ya€[u,t+u]' 

[Y  €[u+u,t+u-hc]'J 


]du 


(t-0) 


x   exp{y 


[-1+P 


u=— °° 


'Y   €[u,t+u]' 
a 

[Yb€[t+u,t+u+x]'J 


]du) 


(t=0) 


(6.7) 


The  second  term  in  (6.7)  becomes  exp(-yx)  when  t  =  0  since  the  exponent 
is  then  of  the  form  (5.5);  the  first  term  may  be  alternatively  written  as 


+CO 


[P(u<Y  <t+u)  +  P(t+u<Y  <t+u+x)  -  P 
a  b 


u= 


u<Y  <t+u 
a 

t+u<Y  <t+u+x 
b 


du] 


(t-0)  ' 


which  then  simplifies  to 


1  -  D 


r+oo 


Ti= 


u<Y   <t+u  1 

a 

t+u<Yt  <t+u+x 

c  b  } 


du 


(t=0) 


(6.8) 
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After  writing  the  derivative  in  the  incremental  way  and  examining  the  expression 
at   t  =  0,   it  is  found  that  (6.8)  becomes  P(Y  -Y   [0,x]').   Thus  from  (6.7), 

D   Si 

we  have  the  main  result  of  this  section'  as 

P{X(b)>x}  =  P(Y  -Y  €[0,x]')exp(-yx),  (6.9) 

a  b  a 

and  it  is  clear  that  this  is  a  proper  survivor  function.   One  particular  case 
of  (6.9)  is  when  Y   and  Y   have  exponential  distributions  of  parameters 
X   and  X  ;   then  the  p.d.f.  of  X     has  the  mixed-exponential  form 

f,   N        X  -(y+X,  )x     X, 

r(b),  N   a  .    , ,  v      b    ,    b    -yx  .,    ... 

fa   (x)  =  rrxT  (v+xb)e  +rrx7  e    •  (6-10) 

a  b  a  b 


The  distribution  of  Y  -  Y   appearing  in  (6.9)  also  features  in  other 

aspects  of  the  bivariate  delayed  Poisson  process,  as  mentioned  by  Cox  and  Lewis 

(1972),  and  the  reason  is  fairly  evident.   For  some  purposes,  the  main  Poisson 

process  can  be  eliminated  from  the  model  by  concentrating  on  (say)  just  the 

type-a  events;  a  type-b  event  occurs  at  a  distance  Y  -  Y   from  its  type-a 

event,  and  these  follow  a  Poisson  process  anyway.   Such  reasoning  leads  to  a 

quicker  derivation  of  (6.9),  but  our  aim  here  has  been  to  apply,  for  the  first 

time  it  is  believed,  the  semi-synchronous  definition  (6.1).   It  should  be  noted 

however,  that  the  univariate  process  with  displacements  Y,  -  Y   is  not  iden- 

b    a 

tical  to  our  delayed  Poisson  process  in  all  respects;  for  instance,  in  terms 
of  initial  conditions  and  numbers  of  events  in  the  course  of  being  delayed. 
Some  joint  counting  distributions  for  the  univariate  delayed  Poisson  process 
have  been  obtained  by  Milne  (1970)  when  pursuing  identif iability  questions. 
Finally,  we  note  that  a  joint  distribution  of  interest  is  that  of  X 

3. 

(a) 
together  with   that   of     X        ,      the   interval   from  the   arbitrary   type-a  event   to 

3. 
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the  next  type-a  event.   We  can  define  generally  the  required  joint  distribution 


as 


P{X(a)>x,  X(b)>y} 
a       a 


lim  P{N  (t,t+x)  =  0,  N.  (t,t+y)  =0   N  (0,t)  ^  1>, 

,  *        3  D  3. 

t-J-0 


(6.11) 


or  use  the  univariate  delayed  Poisson  process  with  Y  =  Y  -  Y   as  the  dis- 
placements.  We  follow  the  latter  course.   Taking  our  origin  as  an  event  in  the 
Poisson  process  of  type-a  events,  we  require  that  there  be  no  type-a  events  in 
(0,x]   and  no  type-b  events  in   (0,y],   We  consider  type-a  events  which  have 
been  generated  in   (-T  ,0]   and   (x,T?]   and  then  let  T   and  TL   tend  to 
infinity;  thus 

P(X(a)>x,  X(b)>y)  = 
a       a 


»  -yT  (yTx)' 
lim  I   e 
T  -k.  r=0 

2 


-Tl 


P(Y€[v,v+y]')dv 


-Mx 


v=0 


°°  -yT_  (yT0)s  n 

L  e        ~rr- 


s=0 


~    P(Y€[-v-x,y-v-x]')dv 
"  2  v=0 


P(Y€[0,y]')    (6.12) 


which  may  be  simplified  to 


0 


P(X(a)>x  X(b)>y)  =  exp{-yx-yy+y    P(u<Yv-Y  <u-Hy)du>  x  P(Y,  -Y  €[0,y]f).   (6.13) 
a      a  J       b  a  ba 

u=-x 
This  joint  distribution  has  the  marginal  distributions  already  discussed.   From 
a  statistical  analysis  point  of  view  it  would  be  useful  to  obtain  the  correla- 
tion between  X     and  X   ;   however,  it  appears  from  (6.13)  that  no  simple 
a         a 

general  expression  will  be  available,  although  (6.13)  should  yield  the  result 
with  particular  assumptions  on  the  delay  structure. 
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